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Abstract 



We study the detailed properties of Z2 domain walls in the deconfined high 
temperature phase of the d = 2 + 1 SU{2) gauge theory. These walls are 
studied both by computer simulations of the lattice theory and by one-loop 
perturbative calculations. The latter are carried out both in the continuum 
and on the lattice. We find that leading order perturbation theory reproduces 
the detailed properties of these domain walls remarkably accurately even at 
temperatures where the effective dimensionless expansion parameter, g'^/T, 
is close to unity. The quantities studied include the surface tension, the 
action density profiles, roughening and the electric screening mass. It is only 
for the last quantity that we find an exception to the precocious success of 
perturbation theory. All this shows that, despite the presence of infrared 
divergences at higher orders, high-T perturbation theory can be an accurate 
calculational tool. 
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1 Introduction 



Non-Abelian gauge theories possess many surprising aspects. An example is 
the hnear confining potential present in both the three and four dimensional 
SU{N) theories at low temperature, T. At some finite T there is a phase 
transition and confinement is then lost|jI|, but this is not unexpected be- 
cause simple energy versus entropy arguments tell us that a confining 'fiux 
tube' will condense into the vacuum at some finite value of the temperature. 
These phenomena are non-perturbative and have so far defied analytic, as op- 
posed to numerical, approaches. However, at sufficiently high T, there would 
appear to be an important theoretical simplification. All these theories are 
asymptotically free so that the effective interaction on the relevant energy 
scale, T, should become small at high T and the physics of the gluon plasma 
should become accurately calculable in perturbation theory. Unfortunately 
there are infrared divergences in higher orders of perturbation theory, which 
are associated with the perturbative masslessness of the magnetic gluon. Al- 
though we expect this gluon to acquire a mass through the non-perturbative 
physics of the dimensionally reduced theory, this leaves room for uncertainty 
about how reliable high-T perturbation theory really is. 

However even this naive picture of gauge theories at high T — as a 
weakly interacting plasma of gluons — contains surprises. There turns out 
to be a symmetry associated with the centre of the group, Z{N), which is 
spontaneously broken at high T. Separating two different Z{N) vacua will be 
a domain wall whose properties have been calculated in perturbation theory 
for T — > oo [Q. However the reality of these domain walls is controversial 
for a variety of reasons. Firstly, there are the general doubts about high-T 
perturbation theory that we alluded to above. Secondly, all this is in the 
usual Euclidean formulation of finite temperature field theories and there 
is a question of what if anything these domain walls might correspond to 
in Minkowski space-time. Finally the walls have peculiar thermodynamic 
properties, which become more acute when one includes quarks into the 
theory ||^. It is important to resolve these uncertainties not only because 
of the theoretical interest of these domain walls, but also because related 
structures can be associated with important physical phenomena when one 
considers the Standard Model in the early universe 0. 

In this paper we address the particular problem associated with the un- 
certain status of perturbation theory at high T. We shall do so by calculating 
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the properties of the domain walls both in perturbation theory and by a fully 
non-perturbative Monte Carlo computer simulation. We shall work with the 
SU{2) gauge group because the problems should not be any different for 
larger groups. Moreover we shall work in 2+1 dimensions rather than in the 
more physical case of 3+1 dimensions. The reason is that the computational 
resources needed are much less in the former case, and only there can we 
perform calculations with enough precision and control to be really useful. 
At the same time, the origin of the infrared problems is similar in d = 2 + 1 
and in c? = 3 + 1, and their severity is, if anything, greater in the lower 
dimensional case. 

We shall find in the calculations described below that one-loop pertur- 
bation theory does indeed work remarkably — even precociously — well for 
our Z2 domain walls. By implication, this provides evidence for the general 
applicability of perturbation theory at high T. 

At the same time we emphasise that we make no attempt, in this paper, 
to address the other controversial aspects of domain walls and so do not at- 
tempt to settle the interesting question of their potential role in, for example, 
separating bubbles of different vacua in the early universe. 

We shall now outline the contents of this paper. In Section 2 we give 
a heuristic introduction to domain walls and the thermodynamics of gauge 
theories at high temperatures. This will provide the general background for 
the more detailed and specific calculations of the later sections. We then 
turn to the perturbative calculation, at one-loop, of the properties of domain 
walls. In Section ^ we perform the calculation of the surface tension for 
the continuum theory. The calculation is by a method that easily extends 
to other dimensions and so, as well as obtaining expressions for d = 2 + 1, 
we can compare with previous d = 3 + 1 calculations. In order to compare 
with our later numerical work it is useful to have similar results for a finite 
value of the lattice spacing a. This calculation is carried out in Section ^. In 
Section ^ we calculate the action density in the domain wall, as a function 
of the distance from the centre of the wall, since this is one of the quantities 
we shall later calculate numerically. In Section |6.1j we calculate the effects 
of oscillations of the wall once its length is large — the 'roughening' of the 
wall, which has been neglected in most previous studies. In Section 6^ we 
show how the finite size of the volume, in the direction orthogonal to the 
wall, affects the profile and surface tension of the wall. This is important to 
understand since our numerical work will necessarily be on lattices of a finite 
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size. We then turn, in Section 7, to a description of our simulations and the 
results we obtain thereby. We begin, in Section 7.1, by describing how we 
can simulate domain walls through the use of twisted boundary conditions. 
In Section 7.2 we point out that these domain walls can be viewed as 't Hooft 
disorder loops that have been squashed by the short Euclidean time direction. 
Section 7.3 describes how we calculate the properties of the domain wall and 
the electric screening mass. Section 7.4 lists the large- volume raw 'data' from 
which we will eventually extract physical quantities. Section 7.5 describes in 
detail how we control finite volume effects. This is crucial since the potential 
problems with perturbation theory are infrared ones. Finally, in Section 7.6 
we compare our numerical results with those of perturbation theory. Section 
8 contains our conclusions. 

Our theoretical analysis, in Sections 3 to 6, is performed for the general 
case of SU{N) gauge fields in 2 < d < 4 dimensions. Our numerical results, 
on the other hand, are for the particular case of SU{2) gauge fields in d = 
2 + 1. The preliminary results of this study appeared in the Proceedings of 
the 1994 Lattice Conference A study of the case of SU{3) gauge fields in 
d = 2 + 1 has recently been reported Both the SU{2) and S'?7(3) results 
are in agreement with our theoretical analysis. In addition both show the 
same large deviation from D'Hoker's self-consistent formula for the Debye 
mass. 

2 General considerations 

In 2+1 dimensions the gauge coupling, (7^, is dimensionful and its value sets 
the mass scale for the theory. If we perform a perturbative calculation of a 
quantity in which there is a dominant momentum scale, Q, then the effective 
expansion parameter will clearly be /Q so that the theory rapidly becomes 
free at short distances. So at high T the effective expansion parameter, 
g'^/T, will be small and we can expect that we should be able to apply 
perturbation theory. All this is very similar to the case in 3+1 dimensions. 
There the coupling is dimensionless but this difference is only apparent: the 
scale-invariance is anomalous, the coupling runs and its value only serves 
to set the overall mass scale (dimensional transmutation). The coupling 
becomes small at short distances, so that its value at high temperature, 
g'^(T), should be small enough for us to apply perturbation theory. In other 
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ways the two theories are also similar: numerical simulations |]2^ show that 
the d = 2 + 1 theory has linear confinement, a deconfining temperature, Tc, 
and a glueball spectrum that are similar in many ways to that of the theory 
in 3+1 dimensions. 

So we expect a hot gauge theory to be, to a very good approximation, a 
plasma of free gluons, with interactions given in terms of the small effective 
coupling at the ambient temperature T. These gluons are screened just as 
are photons in a plasma of charged particles. The Debye screening mass, 
rriD, that they acquire grows with T; to lowest order m^D is 0{g^) so we 
expect, on purely dimensional grounds, that rriE, ~ g(T)T in (i = 3 + 1 and 
tud ~ in = 2 + 1. 

However this simple picture is not the whole story. An additional and 
important role is played by the center, Z{N), of the gauge symmetry group 
SU (N). This sub-group has a special status because the gluons, which trans- 
form according to the adjoint representation, are invariant under gauge trans- 
formations that belong to the centre. Since the gluons only feel SU{N) gauge 
transformations modulo any Z{N) transformation, their symmetry group 
is really SU{N)/Z{N). Sources in the fundamental representation, on the 
other hand, are not invariant under transformations in Z{N). Consider such 
a heavy source in the usual Euclidean space-time formulation of the high-T 
field theory where the Euclidean time is periodic with period 1/T. As is well 
known the presence of such a static heavy source leads to an extra factor of 
P = (l/N) tr Pexp{i / Aodt} — a Polyakov loop — in the partition func- 
tion. In the low T confining phase (P) = while at high T — 1. Since 
the gluons are screened, the correlations are short range at high T and so 
if, for example, P is close to 1 at one point, the vacuum will have P ~ 1 
everywhere. But the Z{N) symmetry tells us that there must be N such 
vacua, in each of which the physics is identical, and which are differentiated 
by P being close to one of the complex A^'th roots of unity. The picking 
of one of these vacua corresponds to the spontaneous breaking of the sym- 
metry. That the spontaneous breaking takes place at high T rather than at 
low T is remarkable but not impossible; it is, for example, a commonplace 
in (self-)dual theories. 

As soon as we have spontaneous symmetry breaking we have the possi- 
bility of domain walls which will occur at the interface between two of the 
vacua. (Throughout this paper we shall call these objects 'walls', and will 
speak of their 'surface' tension, even though the interface is really a string 
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when we are in 2 spatial dimensions. This is to avoid confusion with the con- 
fining string.) In the case of interest to us, SU{2), such a domain wall would 
separate regions of our Euclidean space-time volume that are characterised 
by P ~ +1 on the one side and P ~ —1 on the other. One can compute 
the free energy density of the wall for small coupling and one indeed finds a 
positive excess over the free energy of the gluon plasma. 

Thus we have a picture of the gluon plasma that parallels that of a 
ferro-magnetic substance below the Curie temperature: the average of the 
Polyakov loop, which arises from the heavy fundamental source, being the 
order parameter. But what is the analogue of the external field? Such a 
field is needed in order to make the system choose a particular direction of 
magnetisation on a macroscopic scale. Without such an analogue one cannot 
trigger a thermodynamical state where our order parameter, P, takes the 
value -|-1 or one where it takes the value —1. 

As we remarked in the Introduction, the 'reality' of these domain walls 
is controversial. Of course, the walls have been observed p in = 3 + 1) 
Monte Carlo simulations of the theory, but this has only been for large lattice 
spacings, and there have been speculations that the walls would not survive 
into the continuum limit 0. In this paper we shall address the question of 
whether these domain walls do indeed exist in the continuum limit and, more 
specifically, whether high-T perturbation theory is reliable. But first, in this 
section, we introduce the basic framework within which we work. 

2.1 Yang-Mills fields: the basic parameters 

We consider a pure SU{N) gauge theory in 2 or 3 space dimensions. We 
will concentrate in this subsection on the thermodynamic quantities that 
one can define by enclosing the system in a box of size LyLz- (For the three- 
dimensional case we add an x-direction.) The Hamiltonian for this system 
reads: 

where iE is the canonical momentum for A and Bmn = c^m^n ~ dnA^ + 
i[Am, An]. We use the standard notation for the fields a.s N x N matrices in 
the Lie algebra of the defining representation of SU{N). 

We introduce the free energy F of the system in a heat bath at temper- 
ature T through the Gibbs trace over the physical states of the system. By 
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'physical' we mean that the states obey the Gauss constraint 

yE + i[A,E])\^) = 0. (2) 
The free energy is defined by 

exp (-^) = Trphys exp (^-J^^ ■ (3) 

As is well known, the Gibbs trace can be related to the Feynman path integral 
by 



Trphy.exp (^-^^ = J VAoVAexpi^-^S{A)^ 



(4) 



Here the integration over Aq implements the Gauss law, the action S{A) is 
equal to (1/4) / dxdt trF^^^F^j^, and the temperature enters the formalism 
through the fact that we make the potentials periodic in the Euclidean time 
t with period 1/T. The relationship is fundamental in that it allows us 
to calculate the free energy using the whole panoply of methods available for 
calculating path integrals; such as perturbation theory(see Sections 3 and 4) 
and Monte Carlo methods (see Section 7). 

Throughout this paper we suppose the potentials to be periodic in the 
spatial directions. However, this does not mean that any gauge transforma- 
tion A has to be periodic, only that it should be periodic modulo an element, 
Zk = exp{ik^}, of the centre Z{N) of SU{N). That this is so is easily seen 
from the transformation properties of the potential: 

Af: = A^A.A - tA^d.A. (5) 

If say A{Lz) = A{0)zk, then the Zk will commute with all the matrices in 
(^ and so will disappear from the right hand side, leaving the transformed 
potential still periodic We shall denote such gauge transformations by 
A^. If k = {ky, kz) then this is a gauge transformation which is periodic up 
to Zky in the y-direction and up to z^^ in the 2;-direction. 

The interest of these extended gauge transformations lies in two facts: 

• they leave the Hamiltonian (|l]) invariant; 

• they serve to distinguish subspaces, in the space of physical states, 
which possess a given number, Cy and e^, of electric fiuxes in the y and 
z directions respectively. Clearly this distinction is modulo N . 
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The above notion of electric flux is developed in detail in |T0|. The following 
remarks represent no more than a heuristic outline. Suppose first that we 
have opposite fundamental sources at Xi and X2. To make this system gauge 
invariant we need to join the sources by a finite string P exp{i Adx} 
running between them. Such a string operator creates the unit fundamental 
(electric) flux that must flow between the sources. Suppose we are now in 
the purely gluonic system with no sources and suppose we wish to add a unit 
of electric flux across the whole volume, running in the ^-direction. From 
the above we expect that we can do so by applying to our state the periodic 
string operator tr Pexp{z J^^ dzA^}. Unlike a contractible string operator 
loop, which would represent some local excitation, this operator will clearly 
feel the centre element of the gauge transformation if kz 7^ 0. Indeed it 
is easy to see that it will acquire a factor of Zk^- If we create a state with 
n such units of electric flux in the z-direction, then it will acquire a factor 
of {zkj^. So if we label a state with electric flux e = {ey,ez) by |e), it 
will be an eigenstate of with eigenvalue exp{ik-e^}, (assuming trivial 
transformation properties for the fluxless state |0)). Clearly the state with 
+ 1 fluxes has the same transformation property as that with 1 unit of 
flux; as one would expect in a non-Abelian theory. 

Since the Hamiltonian is invariant under A^, the energy eigenstates can 
be simultaneously labeled by e and we can define a free energy Fe by re- 
stricting the Gibbs trace in (|^) to a given electric flux sector: 

exp = Tre exp (^-J^^ (6) 

With this definition relation (^ can be rewritten as [|10 



exp 



^J^expj-zfc.e^jz^, (7) 



T ) m 



where 



k 



Zk = Trphys exp [ j . (8) 



The Gibbs traces in (||) can be expressed as "twisted" path integrals and 
can be computed using Monte Carlo methods. These path integrals have the 
following defining property. Suppose, for example, that k = (0, fc^). Then 
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due to the occurrence of the transformation Aj^ in (j^) one picks up a factor 
of ex.]i{ikz2Ti /N} in the gauge transform relating A{t, y, 0) to A(t, y, Lz) and 
A[Q,y,z) to A{Lt,y,z), after going around the boundary of the box in the 
z-t directions. This multivaluedness does not affect the gluon field A^. In 
Section 7 we will perform Monte Carlo simulations of such a twisted partition 
function. 



2.2 The two important parameters: string and surface 
tension 

Now we are ready to give a thermodynamic characterisation of the various 
phases of the gauge theory, through the behaviour of the flux free energies 
Fe as the temperature is varied. To obtain simple formulas we will restrict 
ourselves to SU{2), but the end results will be valid for any SU{N). The 
interesting quantity is the flux free energy Fqi in the elongated direction. We 
want to compare it to the flux free energy Fqo and see how the difference 
behaves for low and high temperatures. This is straightforward using for 
the flux free energies. But we still need some theoretical input on how the 
twisted functionals behave. What one finds for low T is that 

l-|^ = Cexp(-^L.), (9) 

while 1 — (Z^i/Zqo) is exponentially smaller for any ky ^ 0. On the other 
hand for high enough T one finds: 



^01 




Dexp )^Ly\Zoo (10) 



and Zkyi is exponentially smaller for any ky ^ 0. The C and D are some 
pre-exponential factors. The evidence for (P) and (10) comes from Monte 



Carlo simulations, as well as analytic Hamiltonian analyses of gauge Potts 



models |]Tl|. So one gets for the free energy difference: 

Foi - Foo ~ p(T)L, if T is small (11) 

Foi - Foo ~ exp \-^Ly \ if T is large (12) 
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Clearly what ( |Tl| ) is telling us is that at low temperatures we are in a confining 
regime, where imposing unit electric fiux across the lattice costs us an energy 
that is proportional to the length traversed by the fiux; and the tension of 
this fiux 'string' is p. So the free energy difference becomes very large at large 
but is insensitive to the transverse spatial dimensions. At some critical 
temperature T^. this behaviour changes into that of (0). The free energy 
difference now becomes exponentially small with the transverse size. This 
behaviour suggests that there is a wall, with an energy density independent of 
the transverse direction, and a total energy proportional to a. This quantity 
has been computed in a semiclassical approximation p, at very high 
temperatures, where perturbation theory should apply. The purpose of the 
present study is to check the validity of this approximation using Monte Carlo 
methods. 

The way the surface tension enters the free energy difference is through 



the exponential. This was first noticed by Bhattacharya et al. |T^ and was 
recently discussed in a Z{2) gauge model [Q. It is reminiscent of an energy 
difference induced by tunneling. As we will see in the next section this is 
indeed a tunneling through a potential that arise from quantum one-loop 
effects. 



2.3 Effective Action and Polyakov Loops 

As we have seen above, the fundamental quantity of interest is the ratio of 
twisted path integrals Z^. These are computed by converting from the vector 
potentials, as integration variables, to Polyakov loops: 



Vl{x) = Pexp<z j dtAQit,x 



(13) 



and by integrating out the remaining variables to get an effective action for 
fl{x). In a suggestive notation: 



exp{-S'efr(ll)} = J'^^^ ^- Pexp ji^^^^dtAo 



exp < — :jS{A) 



9^ 



(14) 

The reader, in looking at this equation, should keep in mind that only the 
eigenvalues Aj of the loop Q are gauge invariant. So it is only these that 
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should appear in the delta function constraint and in Ses- Note also the 
relation between Q and P: 



^'■■"^IfEA.. (15) 



The 5'efr has been worked out |]T3|, |T2], to two loop order in the d = 3 + 1 
case. In Section]^ we will derive the one loop result for any d, and in particular 
for the case of interest in this paper, d = 2 + 1. 

It is important to note that the effective action does not depend on the 
boundary conditions if the volume is large enough. It is also easy to see that 
the twist k = (0, k^) of the previous section corresponds to the following 
boundary conditions in the effective theory for Q{y,z): 

niy,L,)=exp{th^}niy,0). (16) 

With such boundary conditions Zf^ is given by: 

Zj^ = J^^VQexp{-S,sm. (17) 

We also note that the path integral has a formal resemblance to that 
of a spin model partition function. The effective action, when evaluated in 
perturbation theory, will start with a classical kinetic term. The first non- 
zero contribution for constant Q{x) appears at one-loop. So to this order we 
can write: 

' T 



^efr(^^) = / tr 
Jx 



At high T the coefficient of the gradient term is large and we can expect that 
the path integral will be saturated by smooth configurations of Polyakov 
loops. We shall see in the next section that Kg reaches its minima when 
all eigenvalues of fl coincide. Due to a condition det $1 = 1 there are N 
such minima: Q G Z{N). The order parameter distinguishing between these 
degenerate phases is the value of the Polyakov loop P ([T5|). The Z{N) 
symmetry in the effective theory of Polyakov loops is due to the existence of 
gauge transformations which are periodic in the Euclidean t direction up to 
an element of Z{N): A(l/T) = A{0)zk- These transformations leave S{A) 
invariant but multiply Polyakov loops by Zk- 
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In spin model language, we are in an ordered phase when T is large. 
This is a phase in which the entropy of the spin system is low. However 
the spin system is, at the same time, supposed to be describing a very hot 
gauge system with a large entropy in terms of quantum states. It is this 
complementary feature of the spin and gauge systems, that has given rise 
to a lot of confusion in the past few years. For example, one can ask the 
following question: what is the meaning of a localised surface in the spin 
model in terms of the gauge model? In this paper we will not go into this 
question, but take the pragmatic point of view that we just want to calculate 
the exponent in the decay law (JT^) for the hot fluxes. Nonetheless, one 
should bear in mind, that if one does assume the existence of a well localised 
interface, this leads to unusual thermodynamic properties; as we shall see 
below. 

2.4 Some Thermodynamic Properties of the Wall in 
Gauge Theory 

The study of surface effects requires a careful specification of boundary con- 
ditions. Only when that is done can we separate the free energy of the system 
into well defined bulk and surface free energies. Let us suppose we have in 
our two dimensional 'box' a domain wall, separating two domains. Then the 
free energy F will be for large Ly and L^. 

F = fLyL, + aLy (19) 

Here we assume that the size of the box is much larger than any of the 
microscopic quantities in the system. 

What this equation describes is simply a wall with a constant free energy 
density in the y- direction, and a non-trivial free energy profile in the z- 
direction. When integrated over z, this profile gives the interface tension a. 
This profile is like a soliton. 

In our three dimensional SU{2) gauge theory we have a dimensionful 
coupling constant g, with dimension -^mass. The only other dimensionful 
quantity in our problem is the temperature T. We will typically work in the 
regime where g'^/T is a small number. So on the basis of dimensions alone 
the interface tension has a simple form: 

a = a{g^/T)T' (20) 
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with function a positive and dimensionless. We can apply semi-classical 
methods to calculate a. These methods typically give, when applied to soli- 
tons such as monopoles and sphalerons, 

„ scale ^ . V 

^sol 7— 21 

couplmg 

for the energy ii^soi in terms of the scale (Higgs expectation value). So it is 
not surprising that we obtain, as described in detail in the next section, a 
similar result for our profile: 

= "0-5- = "0 • (22) 

7f 9 



where a is a numerical factor of geometrical origin. The expression (|2^) can 
be easily understood if one realizes that the domain wall has a width w of the 
order of the screening length w ~ l/{gT^^'^), carries an excess of free energy 
density A/ ~ and that a ~ wAf. 

We are now ready for the following thermodynamical observations. First 
consider a droplet of radius R of 'minus' phase in a sea of 'plus' phase. The 
free energy excess due to the presence of this droplet equals 

AF = 2TiRa. (23) 

Now, the probability for the appearance of such a droplet is exp{— AF/T}. 
From the explicit dependence of a on T we learn that this probability becomes 
exponentially small for large T. So we find that an ordered phase prevails 
at high temperature^. 

The second observation is a simple consequence of the positivity and 
the explicit temperature dependence of the interface tension. The interface 
entropy equals 

entropy = -—a = '2^0—^ (24) 

So the interface entropy is negative! As is the interface internal energy e, 
although the free energy given by the difference — Tcr is positive. Such a 

study of discrete gauge models reveals the same phenomenon: they can be mapped 
onto a two dimensional Ising model with a coupling ^ logT for T large. That is, for T 
large both models start to order. 
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negative interface entropy is to be expected quite generally for models that 
order at high temperature: the free energy starts to grow at the critical 
temperature, and consequently the entropy is negative. So the presence of 
the wall diminishes the entropy of the system. 



3 Continuum calculation 

In this section we calculate ^efr, the surface tension and the profile of the 
order parameter (Polyakov loop) of the Z{N) wall. When the temperature is 
high the effective gauge coupling becomes small and one can use perturbation 
theory to calculate the effective action (|T^). This has been done by several 
authors ||16|, |1^ for d = 4. Here we derive similar results for d = 3. Our 
method is simpler and is trivially generalized to arbitrary d. The lattice 
version of these results is presented in Section |^. 

We start from the partition function of the pure gauge theory at finite 
temperature: 



dt MrF^,F^, 



= JvAoVAexpl-^l^'^' dtl^tT[{doA + DAoY + B'']^ (25) 



nl/T 

d 

-l/T 



We use some obvious short hand notations. The gauge potentials: = 
iA'J^T"' and the hermitian generators of the SU{N) group are normalized as: 
ti T°'T^ = 6°'''. The covariant derivative is: 

DAo = VAo + t[A,Ao]. (26) 

Bold symbols denote vectors (d — 1 components). The exception is B which 
is an antisymmetric tensor {{d — l){d — 2)/2 components): 

Bik = d,Ak - dkAi + i[Ai, Ak] (27) 



In contrast to [0, we do not begin by fixing the gauge. First we rewrite 
the partition function in the form of an integral over physical fluctuations 
at very high T. These are the transverse components of the vector potential 
A. To achieve this we integrate over Aq which is clearly an auxiliary field. 
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For every given configuration A{x, t) we can do this easily: the integral is 
Gaussian. We obtain: 



X exp |--^ f'^dt I tr [{d^A - DD~^DdoAf + B^] \ (28) 

Thus far our manipulations have been exact. Now we use the smallness 
oi g. In the leading order in g (saddle point approximation) we need to keep 
only terms quadratic in A in the exponent and can neglect the pre-exponent. 
As expected the effect of Aq is to project out the longitudinal fluctuations 
of A from the kinetic term in (pHj). The term does not contain them 
already. We see that the integral over A factorises into integrals over d — 2 
(times N"^ — 1) transversal fluctuations and a trivial infinite factor. 

Each of these integrals represents a very well known partition function of 
the photon gas, the logarithm of which is: 

InZo = {d-2) lndet(-9o^ - V^)-^'^ 
r Vd^~^h 

where V is the volume of space. Thus we obtain the Stefan-Boltzmann law for 
the free energy density of the hot gluon gas (with a multiplicity of 'photons' 
- 1): 

^ = _(iV2_i)^lnZo = -cT'^, (30) 

where c is: 

c=iN'-l)id-2)^^ad). (31) 

Next, we want to find the dependence of the free energy density on the 
value of the Polyakov loop. To this end we calculate the partition function 
( p5D about a constant background field Aq. We make a shift: Aq = A'q + 
and integrate over Aq. This results in simply replacing OqA with 

doA + t[AQ,A] (32) 

in equation (p8|). For a given matrix Aq there will be at least — 1 gener- 
ators that commute with it ('neutral') and at most N{N — 1) that do not 
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('charged'). The 'charges' are given by: Qi — qj, where qi are the eigenvalues of 
the matrix Aq/ {2ttT) . These are related to the eigenvalues A j of the Polyakov 
loop VL ( [T^ ) by: Aj = exp{i27rgj}. The contribution of 'neutral' bosons is un- 
changed by (|32D and is still given by (]29|). The partition function of the 
'charged' gluons is given by: 

In Zq = {d- 2) lndet(-((9o + i2'KqTf - V^^'"^ 

= -{d-2)J ln(l - exp{-|fc|/T + ^2nq}). (33) 

where q = qi — qj- The periodicity in q reflects the Z{N) symmetry. 

Following [|1^ we choose qi = ... = qN-i = q/N. Only — 1 eigenvalues 
are independent: J2iqi = ^ and thus q^ = q/N — q. As q varies from to 1 
the Polyakov loop: 

P=ltrexp|^Ao} (34) 

changes from 1 to exp(i27r/A^). In the SU{2) case q parameterises the only 
path between the two minima of the effective potential for the Polyakov loop. 
In SU{3) it parametrises the lowest action path fl^. It is possible that it 
does the same for > 3. 

With our choice of Aq there are {N — 1)^ 'neutral' gluons and 2{N — 1) 
gluons with 'charges' equal to ±q. The free energy density as a function of 
the Polyakov loop order parameter (parameterised by q) is now given by: 

= _cT'i + V{q)T'^ (35) 

All the dependence on q comes from the 'charged' gluons and is given by the 
universal function: 

V{q) = -2{N -l){d-2)J j^^, ln(l - exp{-|fc|/T}(l - cos 27rg)) 

Vidl2'\ °° 1 

= 2{N -l){d- 2) _(i _ cos27rgz.). (36) 

For d = A this function is expressed via the Bernoulli polynomial B^lq). For 
d = 3 only a numerical evaluation is possible. 



16 



For q varying slowly on the scale of 1/T one can use the tree expression 
for the gradient term [l^. Then the free energy density reads: 



\vqr + Viq)-c 



1 x2 



2 

(37) 

where we defined a dimensionful parameter 7: 



N-1 87r2 

This completes the calculation of Scs of Section 2.3: ^cfr = j^J-' /T. 

The profile of the wall between phases with q{—oo) = and q{oo) = 1 
is given by the function qQ{z) which minimizes (^7p under corresponding 
boundary conditions. It satisfies: 



= M^). (39) 

The width of the wall is controlled by 7: w ~ 7 ~ l/{gT'^/^'^). The exact 
solution qo{z) for d = 4 is known For c? = 3 we solve ( ^91) numerically. 
The surface tension a is given by the integral of the excess of the free energy 
density inside the wall and is proportional to the action on the trajectory 



r rl J rpd/2+l 

' J dz2V{qo{z)) = ^T" dqpV{q) = , (40) 



where we defined: 



«o = 47r J dq^). (41) 



For (i = 3 and = 2 we find using (0): ao = 5.104 . 

An interesting property of = 3 is that the second derivative of V{q) 
diverges at g = mod 1 (see (|36|)). In fact, V{q) oc g^ln(l/g) rather than 
q^ for small q. One can see from ( ^9]) that this results in a Gaussian rather 
than exponential fall-off in the tails of the solution. This is related to the 
fact that the Debye mass in the Thomas-Fermi approximation is divergent 
in (i = 3 if the charged particles are massless: the well-known / rffc^/fc^. On 
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the other hand, D'Hoker argued [|T8[ that all infrared divergences in QCD 



3 



are cut off by a Debye mass m^i of order g^jT\n{T/g'^) (see eq. ( |117D ). This 
means that the fast Gaussian fall-off should saturate at the exponential when 
q gets small enough (27rgT <^ mo) ■, i.e., sufficiently far from the center of 
the wall.0 This behavior of the tail does indeed occur in our Monte Carlo 
study (see Section 7.5). 



4 Lattice version of the interface tension cal- 
culation 

To compare numerical lattice data with theory requires a lattice version of 
the calculation of the previous section. 

In this section we develop a one loop expression on a lattice with Lt sites 
in the temperature direction and an infinite number in all space directions. 
Quantities computed in this section are the effective potential, the profile of 
the wall, and the surface tension. In Section ^ we calculate the expectation 
values of electric and magnetic plaquettes. In Section | the reader will find 
estimates of finite size corrections. In the dimensionality we are interested in 
we have to go to higher loops to get a finite Debye mass. We follow to this 
end the ideas of D'Hoker [|l8l, which amount to a very simple prescription. 

Let us first fix our notation. The continuum action of Section ^ becomes 
on the lattice [pjjj : 



^5cont ^ /?^iat = /3 5] Re[l - ^trt/J (42) 

Here every plaquette p is summed over only once; in contrast to the contin- 
uum action where the sum over fi, v included both orders. 
The lattice spacing a is related to the temperature T by 

a^t = \ (43) 

which follows immediately from the fact that the length of the system in the 
fourth direction is 1/T. 

^ This happens, however, beyond the applicability of our formulas for q{z). 



18 



The lattice action should become the continuum action in the limit where 
the dimensionless quantity aA^^ becomes very small: 

Up ~ exp{ta^F^^,}. (44) 

The lattice coupling becomes in that limit: 

or 

''-M-'=^—^ (46) 



2N ' ~ 



This is the relation between lattice and continuum parameters in any 
number of dimensions d. For c? = 4 all reference to dimensionful proper- 
ties drops out. Perturbation theory is defined by taking the dimensionless 
parameter g'^T'^~'^ small, or /3 large at fixed Lj. 

First we need a definition of the surface tension on the lattice. This can 
be done by taking the lattice analogues of twisted and untwisted boxes and 
their corresponding partition functions. This corresponds precisely to the 
way we measure the surface tension on the lattice (see section 7): 

expj-^L'^-U^^^^. (47) 
I T J Z (untwisted) 

Although our aim is to study a continuum theory, all our Monte Carlo 
results come from a lattice with finite lattice spacing a. The value of a has 
a meaning only in comparison with physical length parameters. For our 
system there are two such parameters: the inverse temperature 1/T and the 
Debye length, or the width of the wall w. The width w ~ 1/ {gT'^^'^~^) is large 
compared to 1/T for small g. To the order we are interested in we can neglect 
corrections due to finiteness of the ratio of w over 1/T or a. We choose large 
P so that the ratio w/a is big: 10 — 20, and we neglect corrections due to its 
finiteness.0 

However, to save computer time we choose to keep the ratio of 1/T to 
a relatively small. This ratio is the number of lattice sites in the temporal 
direction: Lf. To compare with our MC data we calculate the quantities 
discussed in Section ^ for any Lf to the leading order in g. 



We study the corrections due to the finiteness of Lz/w in Section 6.2 
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To achieve this we notice, that to the order we are working in we have a 
theory of free gluons interacting only with the background. All path integrals 
are Gaussian and factorise into integrals over momentum modes. We need 
only to substitute the lattice periodic momenta for continuum ones: 



2 aki J 

— sm = kj : 

a 2 " 

2 . fco + 27rgT ^ 
■sma = fcol^j- 



(48) 



(49) 



a 2 

The lattice momentum varies inside the Brillouin zone: —Ti/a < ki < ix/a. 
The component /cq is discrete and takes values: 27rTn, n = 0, 1, . . . , — 1. 

The partition function for a 'charged' gluon normalized by the 'neutral' 
one is given by: 



^0 



d-2 



lndet(-A;o(g)^ - - lndet(-A;^ - Y^^kj 



(50) 



We use the proper time trick to calculate (0). To simplify formulas we work 
in lattice units a = 1. The r.h.s. of f^Ol) becomes: 



d-2V 
2 T 



fco 



T ^ dki r°° dt 



exp{-(fc2 + EA')0 



(51) 



To rewrite the sum over ko = 0,27r/ Lt, ... ,2tt{1 — 1/Lt) (remember that 
T = 1/L t in our units) we use the Poisson summation formula, which for a 
function /(/cq) periodic with a period 27r has the form: 



E /(^o)= E / ^e^"'^'V(A:o). 



27r(l--l/L0 

E 

feo=0 

We substitute (^2]) into ([51|), shift the variable /cq and obtain for \nZq/Zo: 



(52) 



c^-2 V 
2 T 



E n 



27r Jo t 



exp{-(A;o + iLtuko - i27iqu + J2iki)t} 

- exp{-(A;2 + zLtuko + Eik^)}] ■ (53) 
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Integration over A;^ produces modified Bessel functions: 

— exp{2tcosA; + m/c}, (54) 

and t can be rescaled after tliat. Finally, we obtain the expression which 
replaces V{q), eq. (|36D, at finite Lf. 

fOO rlf <^ 

Vut{q) = 2{N-l){d-2)L'l - ^ e-'^%i,(t)[/o(t)]'-^(l-cos27rgz/) (55) 

It is instructive to see how the continuum limit Lt ^ oo (i.e., aT 0) is 
recovered. For large n and large t>,n the Bessel function behaves as: 

^''Ut)^^=exp{--} and e-%(t) ^ -=. (56) 

Only large t of order {vLtY will contribute to the integral in (|55D. After 
rescaling t by {i^LtY and integrating we get the same expression as in (pG)). 

The variation of q with z can be also taken into account. The gradient 
term on the lattice becomes: 

^{q{z+l)-q{z)r, (57) 



so that the profile equation (|39|) is replaced with: 



g(z + l)=g(z) + ^-^Mat(g), (58) 
where similarly to the continuum case 

f = 8^^^ = 4vr^^^/?^'"'- (59) 

The interface tension is given by an expression similar to (^) except for ao 
being replaced by its lattice version: 



« = 4vrJ^^ J^'dq^VUq) ■ (60) 
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5 Action density profiles 



An important quantity that we measure on the lattice is the expectation 
value of the plaquette action. In the continuum limit this corresponds to 
{E"^ + B^), where the Euclidean E is: 

E = -doA - DAq. (61) 



We also measure separately the expectation values of the plaquettes of each 
orientation, which correspond to (Ey), {El), {B"^) in our d = 3 case. 

An interesting result that we find in our Monte Carlo study is that these 
expectation values display nontrivial profiles, correlated with the position of 
the domain wall. One expects that at high T a perturbative calculation of 
these quantities is possible. In this section we perform such a calculation.^ 



We start again with the partition function (25) on a constant background 
Aq. At high T the effective interaction is weak and we apply a saddle point 
approximation. We write: 




(62) 



where we linearize E and B: 



E = -doA-i[Ao,A]-VAo; B = V x A. (63) 
We have already calculated the integral ( |62D in Section § 

ZiA) ^ Z{q) = exp jjV{q) - c)} . (64) 

Now, to find the action density profiles we want to calculate things like 
{E"^) and (S^), where the average is understood in terms of the probability 
distribution given by the integrand in (^). These averages depend on A, 
or q. The easiest way to calculate {E"^) is to introduce a parameter, say 
e, in front of this term in the exponent, do the integral and differentiate 
the logarithm of the result over e. We shall do this shortly, but before 
that let us improve a bit on the formula (p^). What we need to include 

''Such profiles were also measured in . Here, we show that at high T one can actually 
calculate them analytically. 
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is the contribution of zero point energies of the modes of the fields A. This 
contribution does not depend on the temperature and is not present in (|6^) , 
but it gives an overwhelmingly dominant contribution to {E'^) and (-B^). 
Indeed, this contribution is of order {1/aY, where a is the lattice spacing 
(the UV cutoff), and is much larger than T'^ - the thermal contribution. To 
see what we are missing consider calculating {E'^ + B"^) from (|64D. For that 
one needs only to differentiate (the logarithm of) the right hand side over 
{l/g^). But the right hand side does not depend on g\ This means that 
the quantity {E"^ + B"^) (which is obviously not zero) gives the sum of zero 
point energies and has no thermal contribution at that order. To make this 
fact explicit imagine rescaling g'^ in (|6^) by a factor 1/e. One can absorb 
this factor by rescaling the fields A by ^/t. This will change the measure by 
a factor (y^){Af'-i)('i-i)A/'^ where U = T'^V/a'^ is the total number of the 
lattice sites. The number (A^^ — l){d — VjM is simply the number of non-zero 
modes, or non-zero eigenvalues of the matrix of the quadratic form in the 
exponent (162|) . Eventually we get instead of (|64t): 

Z{q) = (;(^'-i)('^-i)^exp j-T'^-i JjViq) - c)} . (65) 

Using (^5]) we get for any T: 

:{E^ + B'') = ^ ^ a-^. (66) 



2g^ 



It means that each space component of the vector E or the tensor B has the 
average at T = (due to Euclidean invariance): 

This is the dominant contribution and should be compared to /5(1 — trt/p/A^) 
per plaquette which we measure (in our case N = 2, d = 3). It is indeed 
equal to 1 in lattice units up to a small correction. Part of this correction is 
the thermal effect. 

To calculate the thermal part of the action density, multiply only the 
term E'^ in (^) by a factor e. Then consider the following transformation: 
t' = t/ i/e, A'q = v4o-\/e, T' = Ty/e and g'^ = g"^ / y/e. In terms of new variables 



23 



the integral is the same up to a Jacobian factor and thus: 



IX 



exp|-(v/iT)'^-i|^(V(g)-c)} (68) 



g ^{N^-l)id-lW 



Therefore, 



1 /V^ — 1 (7 — 1 

— {E') = {d- l)^- a-' - - V{q))T\ (69) 

The E"^ is a sum of li — 1 components each contributing equally. Using (|6^) 
and (^) we get also: 

1 (d-l)(d-2)N^ -1 , ^^^^ 

— (B^) = ^ L^— a-'^ + - v^(g))r'^. (70) 

Note that while {l/2g^){E^ + B^) is totally due to the vacuum zero point 
energy, it is (1 / 2g'^) {B'^ — E'^) that contains the thermal energy. The vacuum 
term in this quantity cancels in d = 4, because there is the same number of 
E and B components. 

So far we have neglected the variation of q with z. This can be easily 
corrected for. It contributes only to {E'^) the amount (7^/2)(c?2g)^T'^, which 
is equal to V{q)T'^ due to (^). 

Finally, we write down expressions for the plaquette action densities in 
our case {N = 2, d = 3). We use lattice units: a = 1, in which T = l/Lt, 
and P is given by 



^,{B') = l + (c-V^(g))-^ + 0(l//3); 



= l--{c-V{q))j^ + 0{l/(3); 

1 ,„o, 1 1 



^g,{Ei) = l--{c-3V{q))j^ + 0{l/(3). (71) 

The corrections of order 1/(3 are due to non-quadratic terms and are beyond 
our approximation. However, they should be the same for all plaquettes at 
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that order. They cancel in, e.g.: (5^) - {E^) or {B'^){q) - {B^){0). The 
formulas (|7T| ) are in good agreement with our MC data (see Section 7). 

What do we learn from all this? One can see, for example, that the 
thermal part in the fluctuations of (B^) becomes negative at some q near 
1/2, i.e., inside the wall. This means that the {B^) becomes smaller than 
the contribution of the vacuum fluctuations to that quantity !0 This is just 
another side of the old puzzle with negative entropy and thermal energy 
density §1. 

6 Finite size corrections 
6 . 1 Roughening 

Due to the long wavelength thermal fluctuations of its shape, an interface in 
3 spatial dimensions oscillates from its central position by a distance which 
grows as the logarithm of its area. This roughening also occurs in 2 dimen- 
sions, where the effect is proportional to VTy- Here we calculate the effect 
of the roughening on our measurements of the profile of the wall and the 
interface tension and show that this effect is rather small. 

The roughening is due to long wavelength fluctuations of the shape of 
the interface. These fluctuations are therefore essentially classical: ook ^ T. 
We can consider our interface in 2 spatial dimensions as a classical string of 
length Ly with tension and mass per unit length equal to a. The string is 
a set of free oscillators, the normal modes. The amplitudes fk of these are 
the Fourier components of f{y), the shape of the string at a given instance. 
Each oscillator (mode) has energy T in the heat bath (equipartition). On 
the other hand, the mean energy of such an oscillator is cru;^(/|). The mean 
square of the fluctuation of the string is then: 

(/'^c^l/m = (E/.') = E^ (72) 

The dispersion law is: Uk = k. If we take Ly = oo and replace the sum with 
an integral it will diverge linearly: the soft modes get out of hand if Ly does 
not cut them off. For finite Ly the values of momenta q are given by the 

^In other words, (B^) — (Ey) becomes negative. 
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periodic b.c: 

k = —n, n= 1,2,3,..., (73) 
ij 

where k = (translational mode) is removed by our procedure of shifting 
the center of the interface. In principle, our string approximation will break 
down at some large fcmax, but the value of this UV cutoff is not essential for 
the long wavelength effect we are interested in. We put /cmax = oo. Thus we 
have: 

o CO rp T 2 

The factor of 2 is because for each n there is a cosine and a sine mode. The 
sum over n can be evaluated and we get for the mean square deviation of the 
wall from a straight line: 

Aw^ = (fiy)) = (75) 

Let us estimate this effect in our case. Take, for example, Lt = 3 and 
P = 75: 

1/2 




Aw= ^ JLJ50, (76) 



where we used the formula (|40|) for cr and the relation (|45|). For = 12 — 60 
this varies from 0.5 to 1 or so, as compared to w ~ 17. We observe a 
slight variation of the width of the wall of roughly this size in our MC data, 
although it is no doubt optimistic to be applying a string formalism in a 
situation where the length of the wall is comparable to its width! 

We conclude that roughening does not affect our estimates of the profile 
of the wall. This effect is small, because the wall is stiff, and Ly is not too 
large. So what we measure numerically is really the intrinsic profile of the 
interface. 

One can also calculate the correction to the interface tension from the 
string-like fluctuations of the interface using the same idea that the string is 
a set of oscillators in a heat bath. A difficulty lies in the fact that unlike (/^) 
which is a convergent sum of (/|), the amplitudes of the oscillators, the sum 
of their free energies, — Tln(T/c(jfc) + const, is divergent. This divergence is 
ultraviolet, however, and can be subtracted when computing the finite size 
dependence, similarly to the Casimir effect. 
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Another, more illuminating way of deriving this correction is to consider 
the 0-^{q) in (pTl) as an effective potential energy for the long wavelength 
classical thermal fluctuations of the wall. The profile go (-2) satisfying (^Uj) is 
a minimum of J^{q) for the corresponding boundary conditions. The leading 
exponential behavior of the partition function is then: 



^wall 



exp{-^ jjn<lo{x))-m]] = ^w[-^]. (77) 



where we subtracted the bulk free energy. 

The pre-exponential correction to (^) is due to fluctuations of g(?/, z) 
around go (-2) and is given by the determinant: 

det'(-fV2 + \/"(go(^)))-^/^ (78) 

where the prime on 'det' denotes the fact that we omitted the translational 
mode. This mode is proportional to dzQa- Properly normalized it produces 



a factor ^^aLyj^'^T and an integration over the position of the center. The 
spectrum of the operator in (|78|) can be written as 



A = -I'k' + (79) 

where k = 27m/ Ly and A^ are the eigenvalues of — 7^<9^ + V"{qo{z)). The 
A = 7^A;^ band attached to the zero eigenvalue A^ = 0, corresponds to 
fluctuations of the wall as a whole. The roughening effects are due to these 
gapless fluctuations. The corresponding determinant is det'(— 7^^^). We can 
regularize it in the UV by dividing it by a similar determinant with Ly = 00. 
On dimensional grounds: 

det'(-ya;).. Li 

Other eigenvalues are not related to the roughening and we neglect their 
contribution in our estimate. Collecting all the factors we obtain: 



..n - const / dz ^exp {-^} = const ^ ^ -p {-^} . (81) 
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From (^) we can read off the correction to the interface tension^: 



In— f + 0(T/L,) (82) 



2Ly all 



This means that the correction to a, in the way we measure it in our simu- 
lations (see Section 7), is given by: 



Ly y^dp [ T Ly^- ^^^^ 



For example, at /? = 75, Lt = 3 and Ly ~ 50 this correction is only about 
0.2%. 



6.2 Finite corrections 

In the previous section we discussed the roughening effect which introduces 
corrections of the type {In Ly)/Ly to cr. In this section we discuss another 
source of finite size effects: the finiteness of L^. There are two ways the 
finiteness of affects the free energy of the wall. First, there is a correction 
to V{q) itself because it is given by one-loop integrals which depend on 
through the quantization of momenta running in the loop.Q This correction 
should be of order {T /LzY and is relatively small. For L^^w this correction 
is beyond our leading order in g approximation, because {T/wY is of order 

The second source of corrections is due to the finiteness of Lz/w. It 
is obtained by calculating the action of a particle with a 'mass' 7^ which, 
moving in the potential —V{q), returns to its starting value of g in a period 
of 'time' 2Lz This finite size effect we estimate here and show that, for large 
enough Lz, it is exponential in c? > 3, i.e., exp{—Lz/w}, and Gaussian in 
d = 3, i.e., exp{— L^/w^}. 

So we consider a trajectory which starts at rest at q = e and arrives at 
q = 1 — e precisely after a given 'time' Lz so as to satisfy the boundary 

^This correction is analogous to the Luscher's correction ||2^, but for a classical string 
in a thermal bath, rather than a strip (quantum string at T = 0). The universal coefficient 
of the Coulomb correction 1/L to the free energy in Luscher's case becomes in our case 
the coefficient of the In L correction, ft is especially obvious in the interpretation given by 
Stack and Stone 

^There is, of course, a similar dependence on Ly. 
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condition on the Polyakov loop. The trajectory satisfies a Lagrange-Euler 
equation which can be integrated to give ('energy' conservation): 

^^{c^f-V{q) = -V{e)^E. (84) 
We can use it to relate L-^ to e: 

The action can be cast into the form: 

S{e) = £^ dz [^{qr + V{q)] = dq pi\V{q) - V{e)) + V{e)L^{e), 

(86) 

which is convenient for the numerical evaluation of S{e). Also note, that this 
form is familiar in theoretical mechanics as dS = pdq — Edt. The equations 
( |55D and ( PB| ) give a parametric representation of S* as a function of which 
we use to evaluate the correction numerically. 

To find the asymptotic dependence of S on one can use the relation: 

which is a consequence of dS = pdq — Edt, and can be also derived explicitly 
from (H), dH). 

Now we use the asymptotic form of V{q) at small q. For d > 3 it is: 

V{q) = bq' + 0{q') (rf > 3), (88) 



where 6 is a constant which depends on N and d (p6|). Using this form we 
get from (|5B[): 

L,is) = 2^|^\n^ + 0il), (d > 3). (89) 
Integrating the equation: 

^ = be^ + 0{£^) ^ const bexp I- J^lA , (rf > 3), (90) 
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we obtain the asymptotic form of the large correction to S: 



55- - const J^expj-y^L, !>, (d > 3). (91) 



For (i = 3, however, we have: 



where h is given by (]97|). This leads to 



V{q) = hq^\Yi- + 0{q^), {(1 = 3). (92) 



^.(^) = y'^y'ln^ + 0(l), (rf = 3); (93) 

and 

— ^ const ^exp{-^}, (rf = 3); (94) 
Integrating we obtain: 

We see that the asymptotic form of the correction is related to the way the 
tail of the wall decays: exponential in d > 3 and Gaussian in c? = 3.0 
From (|36|) we find for b in (|88D : 

b=i2n)\N~l)id-2)^^ad-2), (rf > 3); (96) 



SS r^- const ^ exp{-— L^}, (rf = 3). (95) 



TC 



while in (i = 3 the value of b defined as in (^) equals: 

6 = 27r(iV-l), {d = 3). (97) 

To get an idea of the size of this correction in our case let us estimate it for 
the case of /3 = 100, Lt = A, L, = 120. We get 7^ = vr^/^Lj = 40071^ ^ 60^. 
The exponent is bLl/A-f'^ ^ (L^/50)^ ^ 5.8 and e"^-^ ^ 1/300. The pre- 
exponent bLz/4 ^ 200. Thus 6S ~ 1. This should be compared to: 

So = dq \f2i^V{^ ^ 50. (98) 

Thus the correction, SS/ Sq, is of the order of a (few) percent. In the following 
section we shall evaluate these corrections, numerically, for all and not 
just large as herein. 

®See, however, footnote |^ on page |l^ and the related discussion. 
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7 Numerical Simulations of Domain Walls 



As we have seen, at high temperatures the theory appears to have degenerate 
vacua which are separated by domain walls. At asymptotic temperatures 
many properties of these domain walls can be calculated in perturbation 
theory; indeed the existence of these interfaces can only be seen when one 
goes beyond tree level. However, as we remarked earlier, the presence of infra- 
red divergences in higher-orders has raised doubts about the applicabihty of 
perturbation theory and, indeed, about the actual existence of the interface. 
To address these doubts we have performed accurate computer simulations 
of the domain walls, and have compared what we find with the results of the 
perturbative calculations. These computer simulations will be described in 
this section. 

If we simulate the high temperature SU{2) gauge theory in a finite but 
large spatial volume, with periodic boundary conditions, then we expect some 
fraction of the field configurations to contain both Z{2) phases in different 
portions of the torus. Such a configuration will contain domain walls sepa- 
rating the two phases and in principle one could study the domain walls by 
focusing on these particular field configurations. However the relative prob- 
ability of such configurations is very small for the temperatures of interest 
and they would not be encountered in a typical Monte Carlo calculation. So 
we have to use an alternative less direct method. What we do is to impose 
twisted boundary conditions on our system, so enforcing the existence of at 
least one domain wall. This will be described in Section 7.1. In Section 7.2 
we show how the domain wall can be interpreted as a 't Hooft disorder loop. 
We then specify the physical quantities that we plan to calculate and describe 
the methods by which we do so in Section 7.3. Section 7.4 summarises our 
Monte Carlo results. Of course, it is crucial to demonstrate that we have 
all finite-volume effects under control — after all, it is infrared effects that 
are the potential problem here — and this we do in Section 7.5. Finally, in 
Section 7.6, we will take our raw 'data' and use it to extract quantities that 
are of direct physical interest in the present context and compare them to 
the perturbative predictions. 
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7.1 Twisted Boundary Conditions 

We work on lattices of size Ly x x Lf in lattice units. The Euclidean 
time extent determines the temperature, aT = 1/Lt, of the field theory. The 
partition function contains the factor exp{—f3S) where f] = 4:/{ag^) and the 
lattice action is as in (|^): 

5 = i:[l-^trf/,], (99) 

where Up is the path ordered product of the SU{2) matrices, Ui, on the links, 
/, that form the boundary of the plaquette p. 

The simplest and most usual way to introduce twisted boundary con- 



ditions is as follows [20|. We change the above action to a twisted action, 
Stw, by replacing trf/p with — trf/p for those plaquettes in the zt plane that 
emanate from the sites {y, z, t) where z and t are fixed to some particular 
values, say z = j and t = k, while y takes all values from 1 to Ly (see Fig. 

1). 

The system with this altered action and with periodic boundary condi- 
tions is equivalent to the system with the original action but with twisted 



boundary conditions |23|. This we see from the following argument. Firstly, 



let us choose labeling z of the sites so that j = L^, i.e. the twist is between 
z = Lz and z = 1. Secondly, to include the possibility of boundary condi- 
tions that are not periodic it is convenient to extend our labeling to include 
2; = 0, as well as z = 1, . . . , L^. If the system is periodic then corresponding 
sites and links with z = and z = are identified (and similarly for other 
directions). Then the system with a twisted action can be viewed as a sys- 
tem with the original action but with fields which are not periodic. To be 
more specific they are periodic except that fory = 1, . . . ,Ly and t = k the 
time-like link at 2; = is mapped into the negative of itself at 2; = L^. This 
is the lattice version of the twisted boundary conditions described in Section 
2.1. 

One can move the zt position of the line of twisted plaquettes by flipping 
the sign of all Ui which bound these plaquettes from one of the sides; but 
one cannot undo the twist completely. It should be clear that the position of 
the twist does not carry any physical significance since it can be moved by 
such a redefinition of the variables Ui. 
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How does the twist lead to the presence of a domain wall? To see this 
consider the same labeling of sites as we have just used. With free boundary 
conditions the system would spend most of the time in one of the two phases 
where Polyakov loops are all near +1 or all near —1. With the twisted 
boundary condition a Polyakov loop at 2; = is mapped onto negative of 
itself a.t z = Lz- Therefore a homogeneous configuration is frustrated and 
the Polyakov loops must create a nontrivial profile in the z direction to 
interpolate between z = and z = L^. 

So to study the high T properties of domain walls we perform Monte 
Carlo calculations on lattices with periodic boundary conditions but with a 
twisted action. How well defined is the domain wall in practice? To answer 
this question we show in Fig. 2a the distribution of Polyakov loops on a 
typical field configuration taken from a 30 x 80 x 2 lattice at /5 = 100. In 
physical units this corresponds to a temperature T ~ 30Tc where is the 
deconfining temperature. We see that the domain wall is very well defined, 
with relatively small fluctuations around a smooth background distribution. 
This is in fact the highest value of T at which we work. The lowest value is 
on a 12 X 30 X 2 lattice at /3 = 7 corresponding to T ~ 2Tc. There a typical 
field configuration looks as in Fig. 2b. The fluctuations are now much larger, 
but the domain wall can still be unambiguously located. So it is clear that, 
for the range of T we study, there is no ambiguity in identifying the domain 
wall. 



7.2 Domain Walls as Disorder Loops 

Before going on to the details of the calculations, we address the following 
natural question. Since the twist is entirely symmetric in z and t why should 
the 'domain wall' separate regions in z rather than regions in t? This question 
can be plausibly answered, in a way that highlights the physics, by first 
considering the twisted system at very low T where the system is manifestly 
rotationally invariant. Here introducing a twist introduces into the system 
a 't Hooft disorder loop |Q which is closed through the boundary in the 



y-direction. This loop will presumably be a flux tube whose width will be on 
the order of the characteristic length scale of the theory, which here is l/g"^. 
Its special property is that it if one considers the gauge potential on a closed 
path that encircles the disorder loop far from its centre, then the presence of 
this loop leads to the potential acquiring a gauge transformation that goes 
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from 1 to a non-trivial element of the centre as we go once around this closed 
path. So if we take a large Wilson loop and pierce it once (or an odd number 
of times) by this disorder loop, then the value of the Wilson loop is changed by 
a factor of —1 as compared to the value it would possess in the absence of the 
disorder loop. Hence the name 'disorder loop'. It is clear that if the vacuum 
contained a condensate of such loops, then these would be sufficient to ensure 
that large Wilson loops varied as the exponential of their area, so that we had 
linear confinement. Suppose we now increase T by reducing the extent of the 
system in the t-direction. Clearly at some point the time extent will become 
smaller than the width of the flux tube, the tube will become squeezed so 
that it extends right across the time direction while still extending over a 
finite region in the 2;-direction. This will occur once T is sufficiently large 
compared to g^] presumably around the deconfining transition. So at high 
temperatures our 'domain wall' is actually a squeezed disorder loop that 
closes upon itself through the y-direction. It is indeed symmetric in z and 
t except for the deformation induced by the limited t extent. The fact that 
the Polyakov loops on either side of the wall have opposite signs is what one 
might expect from such a squeezed 't Hooft disorder loop. 

We have simplified the above argument by assuming that the disorder 
loops exist as definite field fiuctuations in the low-T theory. This is assum- 
ing a great deal of course. Whether they do so exist is one of the central 
questions in the still unresolved problem of colour confinement. This makes 
the connection between these loops and the domain walls at high-T of added 
interest. 

7.3 Quantities Calculated 

We perform calculations on lattices with and without a twist. The simplest 
and most interesting quantity we extract is the extra action, S^, associated 
with the presence of a domain wall. If both twisted and untwisted lattices 
are of the same size, then 

= (5tw) - (^nt) (100) 

where (S't^), ('S'nt) are the average values of the twisted and untwisted actions, 
as defined in Section 7.1. The extra action is related to the free energy of 
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the wall, Fy, = Ftw - -Fnt, by 



|(t)-^» (101) 

where the derivative is taken at constant values of Ly,Lz,Lt. Using this 
relationship we shall test the one-loop prediction for F^. 

It might appear that changing the action for a line of Ly parallel plaque- 
ttes could introduce some additional foca/ contribution to (Stvf) — (Snt) which 
is not related to the free energy of the wall. That this is not so one can see 
by considering a system with 2 parallel twists. This system is equivalent to 
a system without a twist after a redefinition of variables Ui which move the 
twists to a single position where they cancel each other. 

In addition to making predictions for the domain wall free energy, pertur- 
bation theory can also be used to predict the detailed shape of the domain 
wall as it interpolates between the two Z{N) vacua. In the Monte Carlo 
calculations the domain wall is free to move and so if we are to obtain an 
average profile, we need to shift our origin, in each Monte Carlo generated 
configuration, to the centre of the domain wall. We also need to take into 
account the presence of the twist, since the Polyakov loops change sign as 
one moves through it. Our algorithm is as follows. Consider a single Monte 
Carlo generated field configuration. First we average Polyakov loops over y. 
We write this average as p{z). We now want to identify the location, z = Zc, 
of the centre of the domain wall. This is defined operationally as follows. 
We first identify the values of z where p{z) changes from positive to negative 
values (factoring out, of course, the trivial change at the twist itself). Clearly 
the number of such changes must be odd. In practice the domain wall is very 
smooth at high T — as one can see in Fig. 2 — and it is almost always the 
case that there is only one place where there is a sign change. This occurs 
between sites and we shift our origin in z so that the sites where the sign 
changes are labeled hj z = and z = 1. We now ensure that p{z = 1) > 
by multiplying the whole profile by -1, if necessary. Our range of z is now 
from —Lz/2 + 1 to +Lz/2. Somewhere in this range there is the twist and 
the value of p{z) will fiip sign there. If this occurs for z > 1 then we flip the 
signs of p{z) for values of z beyond the twist; if it occurs for 2; < then we 
fiip the signs for z before the twist. In this way we obtain a wall profile with 
p{z > 1) > and p{z < 0) < 0. We can now average this profile over many 
configurations to obtain an average profile. This will be symmetric about 
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^ = 1/2 so we can fold the profile over (with a sign flip) so that it is defined 
for 1 < z < Lz/2 and is positive. This is our final averaged profile. Note 
that in any individual configuration the centre of the domain wall may be 
closer to ^ = than to 2; = 1. That is to say, our profile is 'smeared' over 
distances Sz 1/2. 

In rare cases a given configuration contains more than one sign change in 
p{z) (always factoring out the trivial sign change at the twist). The number of 
these sign changes is clearly odd. There are two possibihties. One is that we 
have a configuration with more than one domain wall, i.e. the one enforced 
by the twist plus pairs that are genuine quantum fluctuations. In this case 
we would typically expect at least one large gap in z between the walls. The 
more trivial possibility is that we might be simply seeing a large fluctuation 
of the values of p{z) near the centre of the wall (where the values are small 
on the average). This would be characterised by very small gaps between the 
locations of the sign changes. The first type of configuration, which we should 
not include in our average, did not occur in any of the calculations that we 
include below. (It does occur if we approach the deconfining transition or 
if we make the extent in y of the lattice, and hence of the wall, sufficiently 
small.) The second type of configuration we should include and we do so by 
taking its centre to be located in the middle sign change. In practice these 
configurations are so rare that there is no visible change in any extracted 
quantities whether we include them or not. 

Having obtained a centre for the wall from the Polyakov loop distribution, 
we can also define an action profile for the wall, and we can clearly do this 
separately for the different fju components of the action. 

A quite different but equally interesting quantity is the electric screening 
mass, tud- This can be obtained from the lightest mass, mp, that couples to 
Polyakov loops, and hence from the tail of the wall profile. We expect that 
for large enough z, 

p(oo) - p(z) oc e-""*^^ (102) 

where p(oo) can be obtained either by working with very large lattices, or 
by performing simulations on a lattice without a twist and using the average 
value of the Polyakov loop obtained therein. So if we define an effective mass 
by 

, , , p(oo) — p(z — 1) , , 

ameff {z) = ]n W \ , 103 
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then 

amp = lim aniesiz) (104) 

In practice we would extract mp from rUes once we were at large enough z 
that the latter had become independent of z. The electric screening mass, 
mpi, should then be given by mp = 2m£i (see below). 

On a finite lattice the above needs to be altered because we expect con- 
tributions going both ways around the 2-torus. So instead of (|102|) we use 

pioo) - p{z) oc e-"*"^" + e-'^"^^(^--^) (105) 

and alter ( |103D correspondingly. 

Since we need to calculate the average action without a twist, we can also 
calculate the screening mass on these untwisted field configurations. Here 
we follow standard techniques for such mass calculations [^]. We construct 
Py = sums of Polyakov loops at each value of z and then obtain the vacuum- 
subtracted correlation function, C{zi — Z2), as a function of their separation, 
z = Zi — Z2. For large separation z we have C{z) oc exp(— amp^;). We define 
an effective mass ames{z) = \ia[C{z — 1)/C{z)], and we increase z until rUes 
becomes independent of z. At this point we can estimate mp = meg. This 
calculation has the advantage that we can prove that meg (-2) > meg (00). In 
practice we modify this formalism for the periodicity in z as described above. 
In addition we calculate with a range of smeared Polyakov loops and use the 
correlation function that minimizes mes{z = 1); in the spirit of a variational 
calculation. However, this turns out not to be really necessary here; unlike 
the situation at T = 0. As we shall see below, this method turns out to be 
much more efficient for the calculation of screening masses than using the 
tails of domain walls. 

7.4 Monte Carlo simulations 

Our Monte Carlo simulations were performed on a variety of periodic lattices 
with and without a twist. We used a standard heat-bath update algorithm 
mixed with over-relaxation steps. 

The control of finite volume effects is particularly important in these cal- 
culations since it is infra-red effects that are usually seen as being at the root 
of any possible breakdown of perturbation theory at high temperatures. We 
have therefore performed extensive numerical checks of finite volume effects 
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and these will be described in detail in Section 7.5. In this section we shall 
confine ourselves to a presentation of those results that have been obtained 
on lattices which are sufficiently large that any finite-volume corrections are 
much smaller than our (very small) statistical errors. This will, of course, 
need to be demonstrated and we shall do so later. 

Now, let us estimate how large the required volumes must be in lattice 
units. If we use a periodic Ly x x Lt lattice, this corresponds to the 
fields in the spatial Ly x volume being at a temperature aT = 1/ Lt. The 
dimensionless inverse coupling /?, as we have seen, is related to the dimen- 
sionful coupling, g'^ by P = 4:/ag^. Thus in physical units the temperature 
is T/(7^ = P/ALt, and, at a fixed value of Lt, f3 (x T. Perturbation theory is 
expected to be most reliable at very high T, so we want to study the theory 
for very large /3. The characteristic length scale at high T is of the order of 
l/gT^^"^: the inverse of the Debye mass l/niD and the width of the wall w 
are of that order. Therefore the spatial sizes in units of a must satisfy: 



From (|106|) we see that for a fixed value of the temperature in physical units, 
T/g"^, the required volumes will be smallest, in lattice units, for Lt = 2 (since 
Lf = 1 is not sensible). Since the perturbative properties of the domain wall 
can be calculated on the lattice, the minimal calculation one might perform 
is to do everything at Lt = 2. However in this case a is as large as possible 
in units of 1/T (aT = 1/Lt = 1/2) and since there have been suggestions 
that high-T perturbation theory might break down as a — we choose to 
perform calculations for several values of Lt. (This will have other advantages 
that will become apparent below.) We shall cover a range of temperatures 
for Lt = 2,3,4 and we shall perform a calculation at one reasonably high 
value of T for the case Lt = 6. In this latter case aT = 1/6 which is surely 
small enough that any breakdown of perturbation theory, as a ^ 0, should 
have become prominent. 

In Table 1 we list the average values of the plaquette, 1 — Snt = ^(trL'p), 
for the calculations without a twist. In Table 2 we do the same for the corre- 
sponding quantity, 1 — St-w, with a twist. We show the values of the lattice 
sizes, the number of Monte Carlo sweeps and the average plaquette action. 
In the twisted case we perform 'measurements' every Monte Carlo sweep; in 




(106) 
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the untwisted case every four sweeps. The typical number of thermahsation 
sweeps prior to taking any measurements is between 25000 and 50000. The 
errors, given in brackets, are typically based on 40 or 50 bins. In a few of the 
lower statistics cases we use as few as 25 bins. The reader will note that at 
some values of the parameters we have several different lattices. These arose 
during the finite volume studies that will be described in detail later on. The 
measurements that we list here are those that do not suffer significant finite- 
size corrections (and are statistically accurate enough to be useful). At the 
different values of Lt we have chosen values of (3 such that the temperatures, 
in units of g"^, are roughly the same, although for higher Lt we are forced to 
cover more limited ranges of T. (The reader may be puzzled that in some 
cases (3 has not been chosen exactly proportional to Lt] for the purposes of 
the work in this paper no particular significance should be attached to these 
choices.) If we now multiply St„ — Snt by the number of plaquettes in the 
twisted lattice, which is the one containing the domain wall, we obtain a 
value for Sy^ and hence, from (|101|) , information on Fw We shall see later 
on what this comparison tells us about the accuracy of perturbation theory. 

As described in the previous section, the lightest mass that couples to 
Polyakov loops, is of particular interest because it is related to the Debye 
screening mass. It can be calculated either from correlations of Polyakov 
loops in the system without a twist, or from the way the tail of the domain 
wall merges into the vacuum once we are far enough away from the centre 
of the wall. In Fig. 3 we show the effective masses as obtained by the two 
methods. In Fig. 3a we have chosen our highest value of T for = 3 while 
in Fig. 3b we show what one obtains for a medium value of T with Lt = 2. 
We see that in both cases the values of amcs{z) as obtained from Polyakov 
loop correlations do become independent of z at larger z, and that these 
'plateaux' occur early enough for the errors to be very small. Since in this 
case mefT(z) is always an upper bound on mp, we can extract an accurate 
estimate of mp using the first value of the effective mass that is, within 
errors, on the plateau. The effective masses obtained from the domain walls 
are clearly consistent with being asymptotic to these mass values. However 
it is equally clear that they would give us much less accurate estimates of 
mp. (We would need to do fits with at least two masses, since there are no 
convenient plateaux, and so the errors on mp would be perhaps an order 
of magnitude greater. Moreover the assumption that the effective masses 
asymptote from below, while reasonable, introduces a difficult to quantify 
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extra systematic error.) So from now on we shall only use the values of mp 
as extracted from Polyakov loop correlations. These are listed in Table 3, for 
those lattice volumes which do not suffer significant finite- volume corrections. 

To obtain the extra action of the domain wall, S^, at a particular value 
of (3, we take the difference St„ — Snt at that /? and multiply by the number of 
plaquettes on the twisted lattice, which contains the domain wall. We expect 
that this extra action will be proportional to the length of the domain wall, 
i.e. to Ly, as long as Ly is not very small. (This and the related question of 
roughening will be addressed when we discuss finite volume corrections.) So 
we form the quantity Sw/Ly, which is the extra action of the wall per unit 
length (in units of the lattice spacing). If we have values of this quantity 
for several values of Ly at a given value of /3 and Lt, we can average them 
to obtain our best overall estimate. In Table 4 we present our final averages 
for this quantity and for the mass anip, as obtained by averaging the values 
given in Tables 1-3. These will form the basic raw material for our later 
comparisons with perturbation theory. 

7.5 Finite Volume Corrections 

We shall be using our values of Sw to test perturbation theory. The details of 
the T and Lt dependence will be important in this comparison. Since the size 
of the domain wall varies with T, it is important that we control any finite 
volume corrections at all values of our parameters. Otherwise part of the T 
dependence we observe might be due to such corrections. In this section we 
describe in detail how we control finite size effects. We begin with effects of 
finite Lz and then consider finite Ly. 

To establish how the finite periodicity in the z-direction affects the action 
of the domain wall, we perform numerical calculations for a large range of 
values of L^. Since these effects may well vary with the lattice spacing, i.e. 
with aT = 1/Lt, we perform such calculations for two different values of Lt, 
but at the same value of physical temperature T/g^. Since the finite-size 
corrections may differ for the contributions that are leading and non-leading 
in jT , we also perform the calculations for two different values of T j at 
the same value of Lt- The parameter values and the corresponding values 
of Syj/Ly are displayed in Table 5. As we discussed previously, see ( p.06| ), 
the natural scale for the domain wall should be of the order of XjagT^I"^ = 
[3Lt/2. This is the scale that appears in perturbation theory (^,|^): 7 = 
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TTy/ f3Lt. We therefore plot in Fig. 4 the values shown in Table 5 against the 
scaled lattice length L2/7. We see that to a good approximation the finite 
size effects are indeed just functions of this scaled length. We also see that the 
finite size effects vary from being very large to being very small over a narrow 
range of values of I/2/7. Indeed, the domain wall effectively disappears for 
L2/7 < 0.8, and while the corrections are still large for ~ 0.9 to 

1.1, they have become invisible, within our statistical errors by the time 
Lz/l ~ 1.35. We therefore see that the values in Tables 2-4, which all 
correspond to lattices satisfying -^2/7 > 1.65, are effectively for = 00. As 
a final precaution against the unexpected, we show in Table 6 some further 
calculations obtained for a wide range of values of (3 and T. Taking these 
together with the values in Tables 1-3 confirms that the pattern of finite size 
effects we see in Fig. 4 is indeed characteristic of the range of T and a covered 
by the calculations in this paper. 

Since the finite size effects appear to be insensitive to the value of T, 
it is interesting to ask whether they can be reproduced in leading order 
perturbation theory, which, after all, is supposedly exact in the T = 00 limit. 
The appropriate formalism is that of the 'ball rolling in the inverse potential' 
as described in Section 6. To that order the finite size effects are functions of 
the scaled length L^/'j, as can be seen from eqs. (|85| , |86D . In Section |0| we 
solved the equations analytically in the limit of large Lz/'~f. In this section 
we solve them numerically for all L^/j, and different Lt, using V[at('?)- We 
note that solution does not exist if the 'time' 2Lz after which the particle 
has to come back oscillating around q = 1/2 is smaller than the period of 
small harmonic oscillations around this point. This gives for the minimal 



value of Lz'. Lz/'y = \J'n/A\n.2 = 1.06 (for Lt = 00). This fits in well with 
what we observe in Fig. 4. For large Lz the correction is exponential in Lz^ 
in this order. This would be difficult to see given our finite statistical errors. 
Moreover we know that in the full theory the correction must ultimately be 
exponential in Lz- From Fig. 4 we see that leading order perturbation theory 
describes the observed finite size effects reasonably well. We also see that 
our criterion, Lz/'~f > 1.65, should be a safe one to use for all values of Lt. 

We now turn to the finite-size effects associated with the transverse length 
of the domain wall, Ly. If Ly is small enough then extra domain walls will be 
produced as quantum fluctuations since the main factor in the suppression 
of domain wall excitation is ~ exp{—j3Sw) and (x Ly. When Ly becomes 
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sufficiently large we expect the leading correction to be that due to roughen- 



ing, as discussed in Section ^]T|. As we see in (^) these corrections should be 
very small for the parameters we use. Of course at asymptotically large Ly 
the profile of the domain wall, defined by averaging the Polyakov loops over 
7/, will broaden as y/Ly. For our values of Ly the broadening of the profile is 



small (see Section 6.1 



To find out what are the corrections at finite L^, we have performed 
calculations for a range of values of Ly. As in our study of the dependence 
we do so for two values of Lt at the same value of T / g"^ and for two values 
of T at the same value of Lt. These are presented in Table 7. We see that 
the surface tension shows no variation with Ly at, say, the 2(T level except 
for Ly = A dX (3 = 25. Here there appears to be a ~ 4 — 8% reduction in the 
tension. This compares well with (^) — recalling that a ~ 6 for = 2, 
although we should certainly not expect (^) to be accurate for such small 

Ly. 

However, while we see that there are no significant Ly corrections to the 
surface tension, this is certainly not the case for the Polyakov loop mass, amp. 
We see that not only does this mass gap show large finite size corrections for 
the smallest values of Ly, but that these corrections become noticeable for 
values of Ly that are not so small. In fact the pattern we see is consistent 
with a relative correction of the form ~ exp(— ampL^). 

The extensive finite-size studies we have carried out in this section show 
that the potentially dangerous infra-red effects are in fact under control and 
that the values of the surface tension and mass gap that we shall be using 
in the next Section, may be regarded as having been obtained on an infinite 
system. 



7.6 Surface Tension and Debye Screening Mass 

Even without looking at the detailed numbers in our Tables, there are two 
properties of the domain walls that are immediately apparent. One is that 
the probability of such a wall being produced at high T is very small. The 
other is that the walls have a finite width. Do these qualitative features 
already teach us something? 

Consider the finite width. This is significant because at tree level a wall 
would have infinite width. That is to say, the width of the wall would be 
oc Lz however large we made L^. One can easily see this by considering the 
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minimum action configuration that interpolates between the vacuum with 
all Polyakov loops +1 and the vacuum with all loops —1. The fact that 
the walls we generate are of a finite width, i.e. independent of L~, once 
is sufficiently large, is implicit in the finite volume studies of the previous 
section. However it is worth showing this explicitly. We define the width of 
the wall, w(90%), as the distance, in lattice units, from the value of z where 
the Polyakov loop is 90% of its asymptotic value to the point where it is of 
the same magnitude but of opposite sign. As a check that there is nothing 
special about the choice of 90%, we shall also define a width, w(2/3), on the 
basis of 2/3 of the asymptotic value. In Fig. 5 we show how w varies with 
Lz in two of the cases where we have made measurements for a wide range 
of lattice lengths. What we see is that the width of the wall does not change 
with increasing once L.^ > 2w(90%): the wall does indeed have a finite 
fixed width. 

How does this width depend on T? We extract w from all our large 
volume calculations (essentially those listed in Table 2) and plot the results 
against 7 = 7r(/3Lt)^/^ in Fig. 6. This variable, as we have seen, determines 
the width of the wall in leading-order perturbation theory. As a matter of 
fact, 7 ~ w(0.97) to that order in the continuum limit, i.e. for Lt = 00. We 
see from Fig. 6 that for each value of Lt the width varies linearly with 7, 
with significant deviations only at the very lowest values of T. The lines for 
different Lt are close to each other, but do not coincide. This is what one 
expects in perturbation theory; the width has to be proportional to the one 
scale, 7 = 2-K/g^/T, but the constant of proportionality will suffer lattice 
spacing corrections, i.e. will depend on aT = 1/Lt. We show in Fig. 6 the 
lines one gets in perturbation theory. Clearly these are consistent. 

The corrections to perturbation theory are governed by g'^/T = 4Lt/f3. 
We note that for /5 = 7 at = 2, and for the corresponding values of (3 at 
other values of Lt, this is > 1. It is therefore remarkable that the calculated 
widths deviate by no more than ~ 10% from the leading-order perturbative 
high T expectations. It would seem that not only does perturbation theory 
work well where we might expect it to, but it even works well where we have 
little reason to hope it might. We shall see other instances of this later on. 

The fact that wall-like quantum fluctuations are very rare tells us that 
the action a wall costs is positive: Sw > 0. Using ( |101| ) and noting that, for 
fixed Lt and Ly, d/df5 ~ d/dT because (3 = A/ag"^ = ALtT/g"^, this tells us 
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that 

A f^U A (h-JIl] > 0. (107) 

dT\T J dT \Lt J - ^ ' 

Here we have used the definition of the domain wall surface tension: 

= aLyCT. (108) 



From (|107|) we immediately deduce that cr{T) oc , where 5 > 2, ignoring 



possible logs. At the same time we expect that it cannot increase faster than 
T^. The perturbative value for the exponent is, as we have seen, 5 = 2.5. So 
we see that our qualitative observations already constrain the temperature 
variation of the surface tension to lie in the interval 6 = 2.5 ± 0.5. In the 
quantitative comparisons below we shall attempt to make the comparison 
with perturbation theory much more precise. 

Before doing so it is interesting to ask what the above bound means 
physically. The following is a simple heuristic interpretation. Let mu be 
the screening mass; then the wall will have a thickness 0{l/mD). A very 
crude expectation is that a oc /nio. Now ii mo grows faster than T then 
it cannot be thermally excited and the whole high-T picture of a screened 
plasma breaks down. The statement that m^) grow no faster than T is 
equivalent to our bound that a grows at least as fast as T^. 

We now turn to a quantitative analysis of the results displayed in Table 
4. We recall that the leading order perturbative result is 

a = -T^-s (109) 
9 

The important energy scale is T and the dimensionless expansion parame- 
ter on this scale is jT . Thus the above leading order perturbative result 
should become exact in the T ^ oo limit. Naively one might expect finite 



temperature corrections to ( 109 ) to be 0((?^/T) - as in 4 dimensions - but 
here in d=2+l there might well be logarithms and the power itself might 
be different. We shall see below that this uncertainty about the functional 
form of the leading corrections will limit the precision with which we can test 
perturbation theory. As we have seen, (|109|) is valid both in the continuum 



and on the lattice, except that in the latter case the value of the constant 
a will receive calculable lattice spacing corrections. Since T is the (largest) 
important physical energy scale, we expect that the lattice spacing correc- 
tions should depend only on aT = 1/Lt. Moreover, since a is a dimensionless 
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physical quantity we expect, on quite general grounds for a pure gauge the- 
ory, that the corrections should be (9(a^T^) ~ Oilj Lt^) for small enough a. 
The detailed perturbative calculations do in fact bear out this expectation. 
From (|101],|108|,|109D we obtain 




to leading order in perturbation theory. We have calculated the value of a 
as a function of aT = l/Lt, using lattice perturbation theory ( |60D and show 
a selection of these values in Table 8. Since the lattice spacing corrections 
depend on Lt, we shall mostly examine the T dependence at fixed Lj. Indeed, 
we have chosen our parameters with this in mind. However before doing so 
we briefly take the alternative approach of varying Lf at fixed (3. For this 
purpose it is useful to rewrite ( |110| ) in the form 

Now, as we see in Table 4, it is only for (3 = 75.0 that we have a usefully large 
range of Lt values. So we plot these values of Sw/Ly against T / = P/4:Lt, 
in Fig. 7. We use logarithmic scales so that a power dependence in T 
will appear as a straight line. And, indeed, the calculated values do fall 
on a straight line to a good approximation. The slope suggests a variation 
oc T^'^ which is close to the perturbative variation of cx T^'^. To carry the 
comparison further we need to take into account the fact that in addition to 
the predicted T^'^ behaviour, there are different aT lattice corrections at the 
different values of Lf. That is to say, the behaviour predicted by perturbation 
theory is oc a{Lt)T^'^ and not just the power of T. In Fig. 7 we show the 
complete leading order perturbative prediction and we see that its variation 
fits that of our data very well. The normalisation is not exactly right — but 
the difference is small and is decreasing with increasing T just as one would 
expect from a higher order correction in g'^/T. If we do indeed try to fit the 
data with a higher order correction we find that an 0{g'^/T) correction will 
not work; but 

a=^T^ ' \ l + 0.13 ( \ (112) 
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fits perfectly well. The precise power of the correction is not to be taken too 
seriously of course; it may be an effective power that partially simulates the 
effects of logarithms in our limited range of T / g"^ (this range being roughly 
3 to 10). 

We conclude from the above comparison that at temperatures T / g"^ > 3 
the T dependence of a is very close to the perturbative expectation of T^-^; 
indeed what we find is ~ T^'^. Moreover this slight difference almost entirely 
disappears when we include the perturbatively calculated aT corrections to 
a. The remnant discrepancy, a few percent, decreases as T increases and so is 
consistent with being a higher order perturbative correction, as, for example, 
in ( p.l2|) . Unfortunately the fact that we do not know the precise functional 
form of this correction, prevents us from carrying out the quantitative com- 
parison any further than this. 

We turn now to consider the bulk of our calculations. We shall consider 
the T dependence at fixed values of Lf = 1/aT, so that the lattice spacing 
corrections do not vary with T. At the same time, by performing calculations 
for several values of aT we can see whether the perturbative predictions show 
any sign of failing as one approaches the continuum limit. To compare our 
results to the perturbative prediction we define a quantity acfr by 



As we see from ( |111D , to leading order in perturbation theory aes(Lt) = 
a{Lt). In Fig. 8 we display our Monte Carlo results for aes as a function 
of g'^/T. The perturbative value of a is also shown, as a horizontal broken 
line, in each case. We observe that the calculated surface tension indeed 
approaches the perturbative value as T increases. At the highest values of T 
the discrepancy is no more than a few percent. Our data is clearly compatible 
with the leading perturbative result being exact in the T — cx) limit. 

In Fig. 9 we plot the ratio aes/a for all our data. We see that to a 
good approximation it is a function only of g'^/T. This supports the idea 
that the small differences we see between the full and perturbative surface 
tensions are in fact due to higher order corrections in g'^/T. We note that 
these higher order corrections are small over our whole range of T. Indeed, 
even for g'^/T ~ 1.1 the correction is only about 25% of the leading term. 
Recall that this temperature corresponds to only about twice the deconfining 




(113) 
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temperature. It is quite extraordinary that lowest order perturbation theory 
should still be so accurate at such low temperatures. 

If we knew the functional form of the leading correction, we would attempt 
to extrapolate our 'measured' values to T = oo. Unfortunately we do not. 
In (i = 3 + 1 we would expect the correction to be simply oc g'^/T. However 
in c? = 2 + 1 there are infrared logarithms which may also resum into a power 
of g'^/T. So it seems reasonable that the correction should be some effective 
power of g'^/T that lies between 0.5 and 1 in our range of T. If we fit our 
data with a form 



then we find that while £ = 1 is excluded, powers near e = 0.65 work perfectly 
well, as we see in Fig. 10. The intercepts of these fits are compatible with 
the leading order perturbative predictions. 

The final question in this context is whether there is any sign that this 
agreement with perturbation theory breaks down as we approach the con- 
tinuum limit. The above detailed comparisons have involved reducing the 
lattice spacing by a factor of 2, i.e. from aT = 1/2 to aT = 1/4. As we 
see in Figs 8,9 there is no sign of any lattice spacing dependence other than 
that calculable in perturbation theory. To go further we also calculated the 
surface tension with aT = l/L^ = 1/6 at = 75. In Fig. 11 we show the 
value of acs/oi{Lt) for this Lf = 6 point as well as for other Lt values at 
approximately the same temperature (using T/g^ = /5/4Lj). Note that since 
the T dependence of this quantity is very weak, as we have seen in Fig. 8, 
it is not important to get T exactly the same at different values of Lt. Note 
also that this is a reasonably high value of T at which to perform such a 
test: the deviations from leading order perturbation theory are only at the 
7% level or so. We see from Fig. 11 that there is no significant deviation 
from perturbation theory with decreasing a even down to aT — 1/6. It there- 
fore seems extremely unlikely that the continuum surface tension will not be 
equally well described by perturbation theory. 

So far we have focused on the surface tension of the domain wall. However 
perturbation theory also makes predictions for more detailed aspects of the 
domain wall, such as its profile. In Fig. 12 we show how these predictions 
compare with our Monte Carlo calculated Polyakov loop profiles for several 
parameter values. We see very good agreement with the main discrepancy 




(114) 
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arising from the fact that the vacuum values of Polyakov loops are ±1 in 
leading order perturbation theory. This is primarily an artifact arising from 
the fact that {9) = does not mean that {cos 6) = 1. A less trivial, but almost 
invisible difference is that the approach to the vacuum at large distances is 
Gaussian for perturbation theory, but exponential in the full theory. 

In Fig. 13 we show some typical examples of profiles of different com- 
ponents of the action density. Again there is very good agreement with 
perturbation theory. 

We turn now to the Debye screening mass, mo- When we expand the 
trace of the Polyakov loop the first non-trivial term is ~ Aq^ so that we expect 
correlations of Polyakov loops to receive contributions from the exchange of 
pairs of screened electric gluons. However, as emphasised by [^|, at higher 



order in the coupling, g'^/T, they also receive contributions from the magnetic 
gluon and larger numbers of gluons of both kinds. The question naturally 
arises: which kind of contribution are we seeing when we extract the masses 
displayed in Table 3? If we expand the normalised correlation function, C{z), 
of Polyakov loop operators, P, in energy eigenstates 

^ = X:cnexp(-E„z) (115) 

then the coefficients, c„, are just the amplitudes squared 

Cn = \{vac\P\n)f (116) 

with normalisation J^^n = 1- From our above discussion we expect to be 
largest for the state with two screened electric gluons. Of course, for large 
enough z, the correlation function will be dominated by the lightest energy, 
Emin, irrespective of the value of Cn (as long as it is non-zero), and so we need 
to check whether the masses we have listed in Table 3 do indeed correspond 
to states for which c„ is large. In Table 9 we list the overlaps, c„, for the 
states whose masses are listed in Table 3. We provide both the overlap onto 
the best smeared Polyakov loop operator and the overlap onto the simple, 
unsmeared Polyakov loop operator. We do so for the extreme values of T 
at each value of Lj. We see that in all cases the normalised matrix element 
squared is > 79%. Given that we expect the magnetic gluon etc contributions 
to receive relative suppressions that are powers of g'^/T, which is ~ 0.1 at 
our highest (3 values, it is clear that the masses we have obtained belong to 
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states that have no such suppression. We therefore claim that we can read 
off rriD from Table 3 using arriD = 0.5amp. 

As we have already seen, the Debye screening mass is infinite at 1-loop, 
due to an infrared divergence. These divergences go away if we use non- 
zero gluon masses in the diagrams and so one can try to do a self-consistent 



calculation for the mass, mo- This has been done by D'Hoker who 
obtains 

In ( —) - 1 + O ^ ^ 



2 9'T 
rriD = 

TT 



.ttidJ ^ ' ^ \iHT/mD))<J\ ^^^^^ 

where ( > 0. We see from ( p.l7| ) that m£,/g'^ is a function only of T/g"^. 
Naively we would, of course, have expected rriD ~ gVT since, as we have 
seen, that is the scale for the domain wall. And indeed this is the leading T 
dependence in ( p.l7| ), up to a weakly varying logarithm. However the correc- 
tion term is down only by logarithms and so we might not be surprised to find 
the comparison with perturbation theory not as good as for the properties 
of the domain wall, where the corrections are powers of g'^/T. 

In Fig. 14 we plot our masses against T/ g"^. What we actually choose to 
plot is niD^ I g'^T since that way we factor out the supposedly dominant g'^T 
factor, and so expose the remaining variation more clearly. We also show 
the leading perturbative prediction as obtained from the first term of 



The first observation is that the dominant variation of m/j is indeed ~ gT^'"^. 
However there is a substantial additional variation which is too strong to be 
due to corrections that are higher order in g^ jT . Indeed if we try a fit of the 
form ^ 

mu' = g'T[c^^ciQ\ \ (118) 



we find that it simply does not work, even if we remove the lowest T point 
from the fit. Indeed, as we see, this additional T variation is quite similar 
to that obtained from ( |1171 ). However the normalisation is completely off. 
Even allowing for the Lt dependence in our values of m^), there is a discrep- 
ancy of about a factor of 3 with perturbation theory. Since the corrections 
in ( |117| ) are only logarithmic we cannot, of course, claim a contradiction 
with perturbation theory. However it is worrying that there is no trend to- 
wards a reduction of the discrepancy even at our highest values of T where 
1/ ln(T/5f^) ~ 0.4. This is in stark contrast to other properties of the domain 
wall where we found the corrections to leading order perturbation theory to 
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be small even for /T ^ 1. 



8 Conclusions 

In this paper we have carried out extensive perturbative and Monte Carlo 
calculations of the high-T domain walls which are associated with the sponta- 
neous breaking of a Z{N) symmetry in SU{N) gauge theories. As we argued, 
these walls can be viewed as 't Hooft disorder loops which become squeezed 
once the extent in Euclidean time becomes small enough, as it does at high 
enough T. Our purpose has been to test high-T perturbation theory and to 
estabhsh whether these unusual objects do really exist in the Euchdean con- 
tinuum theory. In order to be able to obtain numerical results of sufficient 
accuracy to be convincing, we have worked with the simplest theory that 
one may consider as realistic in this context: the SU{2) gauge theory in 2-1-1 
dimensions. 

This kind of calculation is difficult for several reasons. Firstly the po- 
tential problems are infrared and it is therefore crucial to make sure that 
the volumes used are large enough. This means not only doing detailed nu- 
merical finite-size studies, but also calculating the appropriate finite- volume 
corrections in perturbation theory. For example, in this paper we have shown 
how one can calculate the effects of roughening on these domain walls. Sec- 
ondly one-loop perturbation theory becomes exact, at best, only in the limit 
T — s> oo. Now at fixed aT the size of the lattice will grow as T^/^, in lattice 
units, simply because the width of the waU is 0(1/TV2) in = 2 + 1. This 
makes it difficult to simultaneously get close to the continuum limit, where 
aT is small, and to reach very high values of T. For this reason we have 
performed the perturbative calculations not only for the continuum theory 
but also for the lattice theory. In this way we can directly compare pertur- 
bation theory to the full non-perturbative results one gets from Monte Carlo 
simulations. 

In practice we have carried out simulations for temperatures as high as 
~ SOTc, where Tc is the deconfining temperature, and for lattice spacings 
as small as 1/6T. (This comparison with 1/T is appropriate, since T is 
the largest important physical energy scale in the problem). At the highest 
values of T our numerically obtained values of the surface tension agree with 
the perturbative predictions at the percent level and this is so at all our 
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values of a. Moreover there is agreement at the 25% level or so, even at 
temperatures as low as g^/T ~ 1. When we look at the variation with a, 
over the range a = 1/2T to 1/6T, we again find excellent agreement with 
perturbation theory, with not the slightest hint of any anomaly developing 
as a — > 0. 

At the same time we have obtained perturbative predictions for the more 
detailed properties of the wall, such as the action density profile. These cal- 
culations agree very well with our simulations. These profiles are interesting 
in themselves and show, for example, that {B"^ — E"^) — the thermal energy 
(in Euclidean space) — becomes negative inside the wall (see Section |^). 

The only quantity where we fail to find agreement with perturbation 
theory is for the Debye screening mass. However here the perturbative cal- 
culations are not straightforward and the corrections are expected to be a 
power of l/log(T/5f^) and not of g'^/T. So there is no good reason to read 
too much significance into this particular discrepancy. 

Our conclusion is that these high-T domain walls are present in the Eu- 
chdean theory, exactly as predicted by perturbation theory, both on the 
lattice and in the continuum. Since these walls are quantum rather than 
semiclassical objects, they provide a severe testing ground for high-T pertur- 
bation theory. Its success here lends strong support to the usual pragmatic 
assumption that perturbation theory reliably describes gauge theories at high 
temperatures. 
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Table 1: Average plaquette on lattices without a twist. 
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Table 2: Average plaquette on lattices with a twist. 
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Table 3: Masses obtained 



from correlations of Polyakov loops. 
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30x80 


0.2226(42) 


6 


75.00 


36x90 


0.1541(23) 



Table 4: The action density of 
Polyakov loop masses. 



the domain wall per unit length and averaged 



Lt 


13 




a TUp 


2 


100.0 


0.055474(183) 


0.315(7) 




75.0 


0.064096(365) 


0.352(4) 




50.0 


0.08048(22) 


0.4123(25) 




25.0 


0.11773(38) 


0.5269(27) 




15.0 


0.1615(20) 


0.632(8) 




7.0 


0.2598(29) 


0.771(7) 


3 


112.5 


0.02728(18) 


0.2188(44) 




75.00 


0.03342(30) 


0.2572(32) 




37.47 


0.04811(76) 


0.3372(44) 




22.45 


0.06448(119) 


0.395(6) 




10.27 


0.1041(25) 


0.489(4) 


4 


99.97 


0.01788(34) 


0.1897(18) 




49.95 


0.02604(56) 


0.2420(19) 




29.91 


0.03595(80) 


0.2802(23) 




13.81 


0.05661(119) 


0.3411(40) 


5 


62.40 




0.177(7) 




37.33 




0.223(4) 


6 


75.00 


0.01104(41) 


0.1541(23) 



Table 5: The action density of the domain wall as a function of the lattice 
length, Lz, for selected values of (3 and Lf. 



/5 


Lt 




Sw/ Ly 


25 


2 


60 


0.11840(101) 






48 


0.11722(96) 






40 


0.11719(71) 






30 


0.11792(84) 






26 


0.11497(77) 






20 


0.06947(84) 


50 


2 


60 


0.07979(36) 






40 


0.07877(36) 






34 


0.07169(35) 






30 


0.04981(52) 






26 


0.01735(44) 


75 


3 


100 


0.03340(39) 






80 


0.03352(32) 






60 


0.03239(22) 






54 


0.02940(35) 






50 


0.02354(44) 






46 


0.01577(44) 






40 


0.00640(60) 



Table 6: Some additional values of the action density of the domain wall, for 
values of smaller than in Table 2. 









Sw/ Ly 


2 


100.0 


60 


0.05498(18) 




75.0 


50 


0.06391(28) 




15.0 


24 


0.1603(14) 


3 


112.5 


60 


0.02701(16) 




37.47 


40 


0.04863(66) 


4 


99.97 


80 


0.01688(26) 




49.95 


60 


0.02054(58) 




29.91 


48 


0.03557(82) 



Table 7: Variation of the action density of the domain wall and Polyakov 
loop masses with the length, Ly, of the wall. 



/5 


Lt 


Ly 


Sw/ Ly 


a nip 


25 


2 


30 


0.11765(59) 


0.529(3) 






20 


0.11820(86) 


0.529(5) 






12 


0.1184(10) 


0.529(6) 






8 


0.1162(15) 


0.492(11) 






4 


0.1091(21) 


< 0.384(9) 


50 


2 


50 


0.08076(34) 


0.396(12) 






40 




0.416(4) 






30 


0.08068(53) 


0.410(10) 






20 


0.07991(61) 


0.409(6) 






16 


0.08005(58) 


0.386(15) 






12 


0.07937(64) 


0.375(9) 






8 


0.07860(65) 


0.335(8) 






4 


0.07870(179) 


0.272(12) 


75 


3 


120 


0.03271(54) 








60 


0.03409(47) 








30 


0.03347(28) 








24 




0.251(13) 






18 




0.251(7) 






12 


0.03303(65) 


0.218(7) 






6 




< 0.183(6) 



Table 8: Perturbative values of the constant, a, in the interface tension. 



Lt 


a 


2 


6.024 


3 


5.655 


4 


5.409 


5 


5.284 


6 


5.221 


8 


5.165 


10 


5.142 


20 


5.113 


oc 


5.104 



Table 9: Overlaps of Polyakov loop operators onto the states coresponding 
to our values of mr,. 



Li 


3 




Om=i 


2 


100.0 


0.88 


0.81 




25.0 


0.93 


0.82 




7.0 


0.98 


0.90 


3 


112.5 


0.84 


0.79 




10.27 


0.98 


0.89 


4 


99.97 


0.92 


0.81 




13.81 


0.97 


0.86 


6 


75.00 


0.89 


0.80 



Figure 1: Location of the twist. The plaquettes that are indicated will appear 
with a factor of —1 in the twisted action. 




Figure 2: Values of Polyakov loops on typical field configurations with a 
domain wall: for = 2 at: (a) p = 100, (b) /3 = 7. 




Figure 3: Effective masses from Polyakov loop correlations (+) and from the 
tails of domain walls (o) for (a) P = 112.5 with = 3, (b) /? = 25 with 




Figure 4: Dependence of the domain wall action density on the (scaled) 
length of the lattice. Plotted are the values in Table 5 and the lines are the 
leading order perturbation theory expectation (see text). 




Figure 5: The width of the domain wall as a function of the length of the 
lattice for (a) /? = 75 and Lt = 3, (b) /? = 25 and Lt = 2. Widths are 
calculated when the Polyakov loop attains 2/3 (o) and 9/10 (+) of its vacuum 
value. 
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Figure 6: Widths of the domain wall (defined as in Fig. 5) plotted against 
7 = 7r{PLtY^'^ for (a) Lt = 2, (b) Lt = 3, (c) Lt = 4. Lines are corresponding 
perturbative predictions; the dotted lines in (c) are for the continuum, — 
oo. 
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Figure 7: Action density of the domain wall (o) for Lt — 2,3,6 at P — 75. 
The hne is a fit (see text). Also shown are the perturbative values (+) for 
Lt = 2,3,4,5,6. 
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Figure 8: The interface tension (with the factor T'^'^/g removed): numerical 
values (<>) compared to leading order perturbation theory (dashed lines). 



10 



(a) Lt = 2 



to 







1 1 1 1 


1 1 1 






1 1 1 1 


1 1 1 



0.2 0.4 0.6 0.8 1 1.2 1.4 



to! 

CD 




0.2 0.4 0.6 0.8 1 1.2 1.4 



to 

CD 




0.2 0.4 0.6 0.8 1 1.2 1.4 



Figure 9: Ratio of numerical and perturbative interface tensions versus 
high-T expansion parameter, g^/T. 
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Figure 10: Trial extrapolations to the T — oo limit of the calculated interface 
tensions. 
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Figure 11: Ratio, at fixed T, of numerical and (leading-order) perturbative 
interface tensions for decreasing lattice spacing (a = l/L^T). 




Figure 12: Profile of the Polyakov loop within the domain wall (centred at 
z — 0): numerical values (solid lines) and perturbation theory (dashed lines) . 




Figure 13: Vacuum-subtracted local action densities within the domain wall: 
numerical calculations at /5 = 99.97 and Lt = 4, compared with perturbation 
theory (hnes). 




Figure 14: Numerically calculated values of the Debye mass, compared to 
D'Hoker's self-consistent perturbative prediction. 
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